Star Spot Induced Radial Velocity Variability in LkCa 19

We describe a new radial velocity survey of T Tauri stars and present the first results. Our search is motivated by an interest in detecting massive young planets, as well as investigating the origin of the brown dwarf desert. As part of this survey, we discovered large-amplitude, periodic, radial velocity variations in the spectrum of the weak line T Tauri star LkCa 19. Using line bisector analysis and a new simulation of the effect of star spots on the photometric and radial velocity variability of T Tauri stars, we show that our measured radial velocities for LkCa19 are fully consistent with variations caused by the presence of large star spots on this rapidly rotating young star. These results illustrate the level of activity-induced radial velocity noise associated with at least some very young stars. This activity-induced noise will set lower limits on the mass of a companion detectable around LkCa 19, and similarly active young stars.Comment: ApJ accepted, 27 pages, 12 figures, aaste

Observation of enhanced X-ray emission from the CTTS AA Tau during a transit of an accretion funnel

AA Tau was observed for about 5h per XMM orbit (2 days) over 8 successive orbits, which covers two optical eclipse periods (8.2 days). The XMM optical/UV monitor simultaneously provided UV photometry with a ~15 min sampling rate. Some V-band photometry was also obtained from the ground during this period in order to determine the dates of the eclipses. Two X-ray and UV measurements were secured close to the center of the eclipse. The UV flux is the highest just before the eclipse starts and the lowest towards the end of it. We model the UV flux variations with a weekly modulation (inner disk eclipse), plus a daily modulation, which suggests a non-steady accretion. No eclipses are detected in X-rays. For one measurement, the X-ray count rate was nearly 50 times stronger than the minimum observed level, and the plasma temperature reached 60 MK, i.e., a factor of 2-3 higher than in the other observations. This X-ray event, observed close to the center of the optical eclipse, is interpreted as an X-ray flare. We identify the variable column density with the low-density accretion funnel flows blanketing the magnetosphere. The lack of X-ray eclipses indicates that X-ray emitting regions are located at high latitudes. Furthermore, the occurrence of a strong X-ray flare near the center of the optical eclipse suggests that the magnetically active areas are closely associated with the base of the high-density accretion funnel flow. We speculate that the impact of this free falling accretion flow onto the strong magnetic field of the stellar corona may boost the X-ray emission (abridged).Comment: 17 pages and 9 Figures. Accepted by A&

Constraints on the disk geometry of the T Tauri star AA Tau from linear polarimetry

We have simultaneously monitored the photometric and polarimetric variations of the Classical T Tauri star AA Tau during the fall of 2002. We combine these data with previously published polarimetric data covering two earlier epochs. The phase coverage is complete, although not contiguous. AA Tau clearly shows cyclic variations coupled with the rotation of the system. The star-disk system produces a repeatable polarisation curve where the polarisation increases with decreasing brightness. The data fit well with the model put forward by Bouvier et al. (1999) where AA Tau is viewed almost edge-on and its disk is actively dumping material onto the central star via magnetospheric accretion. The inner edge of the disk is deformed by its interaction with the tilted magnetosphere, producing eclipses as it rotates and occults the photosphere periodically. From the shape of the polarisation curve in the QU-Plane we confirm that the accretion disk is seen at a large inclination, almost edge-on, and predict that its position angle is PA~90 deg., i.e., that the disk's major axis is oriented in the East-West direction.Comment: Astron. Astrophys., in pres

Global 3D Simulations of Disc Accretion onto the classical T Tauri Star V2129 Oph

The magnetic field of the classical T Tauri star V2129 Oph can be modeled approximately by superposing slightly tilted dipole and octupole moments, with polar magnetic field strengths of 0.35kG and 1.2kG respectively (Donati et al. 2007). Here we construct a numerical model of V2129 Oph incorporating this result and simulate accretion onto the star. Simulations show that the disk is truncated by the dipole component and matter flows towards the star in two funnel streams. Closer to the star, the flow is redirected by the octupolar component, with some of the matter flowing towards the high-latitude poles, and the rest into the octupolar belts. The shape and position of the spots differ from those in a pure dipole case, where crescent-shaped spots are observed at the intermediate latitudes. Simulations show that if the disk is truncated at the distance of 6.2 R_* which is comparable with the co-rotation radius, 6.8 R_*, then the high-latitude polar spots dominate, but the accretion rate obtained from the simulations is about an order of magnitude lower than the observed one. The accretion rate matches the observed one if the disk is disrupted much closer to the star, at 3.4 R_*. However, the octupolar belt spots strongly dominate. Better match has been obtained in experiments with a dipole field twice as strong. The torque on the star from the disk-magnetosphere interaction is small, and the time-scale of spin evolution, 2 x10^7-10^9 years is longer than the 2x10^6 years age of V2129 Oph. The external magnetic flux of the star is strongly influenced by the disk: the field lines connecting the disk and the star inflate and form magnetic towers above and below the disk. The potential (vacuum) approximation is still valid inside the Alfv\'en (magnetospheric) surface where the magnetic stress dominates over the matter stress.Comment: 15 pages, 15 figures, after major revision, added 3 figures, 2 tables. Accepted to MNRA

A Differentially Private Weighted Empirical Risk Minimization Procedure and its Application to Outcome Weighted Learning

It is commonplace to use data containing personal information to build predictive models in the framework of empirical risk minimization (ERM). While these models can be highly accurate in prediction, results obtained from these models with the use of sensitive data may be susceptible to privacy attacks. Differential privacy (DP) is an appealing framework for addressing such data privacy issues by providing mathematically provable bounds on the privacy loss incurred when releasing information from sensitive data. Previous work has primarily concentrated on applying DP to unweighted ERM. We consider an important generalization to weighted ERM (wERM). In wERM, each individual's contribution to the objective function can be assigned varying weights. In this context, we propose the first differentially private wERM algorithm, backed by a rigorous theoretical proof of its DP guarantees under mild regularity conditions. Extending the existing DP-ERM procedures to wERM paves a path to deriving privacy-preserving learning methods for individualized treatment rules, including the popular outcome weighted learning (OWL). We evaluate the performance of the DP-wERM application to OWL in a simulation study and in a real clinical trial of melatonin for sleep health. All empirical results demonstrate the viability of training OWL models via wERM with DP guarantees while maintaining sufficiently useful model performance. Therefore, we recommend practitioners consider implementing the proposed privacy-preserving OWL procedure in real-world scenarios involving sensitive data.Comment: 24 pages and 2 figures for the main manuscript, 5 pages and 2 figures for the supplementary material

Global 3D Simulations of Disc Accretion onto the classical T Tauri Star BP Tauri

The magnetic field of the classical T Tauri star BP Tau can be approximated as a superposition of dipole and octupole moments with respective strengths of the polar magnetic fields of 1.2 kG and 1.6 kG (Donati et al. 2008). We adopt the measured properties of BP Tau and model the disc accretion onto the star. We observed in simulations that the disc is disrupted by the dipole component and matter flows towards the star in two funnel streams which form two accretion spots below the dipole magnetic poles. The octupolar component becomes dynamically important very close to the star and it redirects the matter flow to higher latitudes. The spots are meridionally elongated and are located at higher latitudes, compared with the pure dipole case, where crescent-shaped, latitudinally elongated spots form at lower latitudes. The position and shape of the spots are in good agreement with observations. The disk-magnetosphere interaction leads to the inflation of the field lines and to the formation of magnetic towers above and below the disk. The magnetic field of BP Tau is close to the potential only near the star, inside the magnetospheric surface, where magnetic stress dominates over the matter stress. A series of simulation runs were performed for different accretion rates. They show that an accretion rate is lower than obtained in many observations, unless the disc is truncated close to the star. The torque acting on the star is about an order of magnitude lower than that which is required for the rotational equilibrium. We suggest that a star could lose most of its angular momentum at earlier stages of its evolution.Comment: 11 pages, 13 figures, submitted to MNRA

On detectability of Zeeman broadening in optical spectra of F- and G-dwarfs

We investigate the detectability of Zeeman broadening in optical Stokes I spectra of slowly rotating sun-like stars. To this end, we apply the LTE spectral line inversion package SPINOR to very-high quality CES data and explore how fit quality depends on the average magnetic field, Bf . One-component (OC) and two-component (TC) models are adopted. In OC models, the entire surface is assumed to be magnetic. Under this assumption, we determine formal 3{\sigma} upper limits on the average magnetic field of 200 G for the Sun, and 150 G for 61 Vir (G6V). Evidence for an average magnetic field of ~ 500 G is found for 59 Vir (G0V), and of ~ 1000 G for HD 68456 (F6V). A distinction between magnetic and non-magnetic regions is made in TC models, while assuming a homogeneous distribution of both components. In our TC inversions of 59 Vir, we investigate three cases: both components have equal temperatures; warm magnetic regions; cool magnetic regions. Our TC model with equal temperatures does not yield significant improvement over OC inversions for 59 Vir. The resulting Bf values are consistent for both. Fit quality is significantly improved, however, by using two components of different temperatures. The inversions for 59 Vir that assume different temperatures for the two components yield results consistent with 0 - 450 G at the formal 3{\sigma} confidence level. We thus find a model dependence of our analysis and demonstrate that the influence of an additional temperature component can dominate over the Zeeman broadening signature, at least in optical data. Previous comparable analyses that neglected effects due to multiple temperature components may be prone to the same ambiguities.Comment: 18 pages, 11 figures, accepted for publication in Astronomy & Astrophysic

Enhanced Perturbative Continuous Unitary Transformations

Unitary transformations are an essential tool for the theoretical understanding of many systems by mapping them to simpler effective models. A systematically controlled variant to perform such a mapping is a perturbative continuous unitary transformation (pCUT) among others. So far, this approach required an equidistant unperturbed spectrum. Here, we pursue two goals: First, we extend its applicability to non-equidistant spectra with the particular focus on an efficient derivation of the differential flow equations, which define the enhanced perturbative continuous unitary transformation (epCUT). Second, we show that the numerical integration of the flow equations yields a robust scheme to extract data from the epCUT. The method is illustrated by the perturbation of the harmonic oscillator with a quartic term and of the two-leg spin ladders in the strong-rung-coupling limit for uniform and alternating rung couplings. The latter case provides an example of perturbation around a non-equidistant spectrum.Comment: 27 pages, 18 figures; separated methodological background from introduction, added perturbed harmonic oscillator for additional illustration, added explicit solution of deepCUT equation